Unifying explanation for carrier relaxation anomaly in gapped systems 
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We develop a theory to describe energy relaxation of photo-excited carriers in low-temperature 
ordered states with band gap opening and formulate carrier relaxation time r near and below 
transition temperature T c by quantifying contributions from different carrier-phonon scatterings 
to the relaxation rate. The theory explains anomalous experimental observations of r in gapped 
systems. Transverse acoustic (TA) phonon modes play a crucial role in carrier relaxation; their heat 
capacity determines r-divergence near T c . The theory is validated by fitting r of fullerene polymers 
onto a theoretical curve. 
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The nonequilibrium dynamics of photo-excited carri- 
ers in solids has attracted considerable research interest 
in the field of condensed matter physics. These dynamics 
are governed by multiple scatterings of carriers and en- 
ergy transfer to the phonon field, both of which are quan- 
tified by the carrier relaxation time r exp and its tempera- 
ture (T) dependence. Usually anomalous T-dependence 
of T cxp is observed in many gapped systems, for example, 
phase-ordered systems showing an energy gap opening in 
the electron band below the transition temperature T c . 
It has been found that a wide variety of superconductors 
[lUlijl and density- wave compounds 13 - show diverg- 
ing behavior of T oxp near T c , as confirmed by femtosecond 
time- resolved optical spectroscopy [H 13. IBl m- fTR Hij . The 
divergence of T cxp in these gapped systems is believed to 
result from recursive energy transfer between electrons 
and phonons. Photo-excited electrons having high energy 
emit a number of phonons through relaxation from above 
to below the energy gap. Conversely, relaxed electrons 
can be re-excited above the gap by absorbing phonon 
energy. This phonon emission-reabsorption process be- 
comes efficient near T c because of the small gap energy, 
and thus, it suppresses the relaxation of carriers, extend- 
ing r exp significantly. This anomalous phonon-mediated 
relaxation in gapped systems is called the phonon bottle- 
neck effect. This bottleneck enables the reproduction of 
these experimental observations of r cxp in various gapped 
systems that exhibit r oxp -divergence at T c . 

There exists, however, a distinct class of gapped sys- 
tems such as Tl-based superconductors 0, Q and Ceo- 
related materials 0, 18 1 that, instead of showing r cxp - 



divergence at T c , show monotonic increases in r oxp with 
cooling across T c . In light of the bottleneck, the mono- 
tonic variation in r exp near T c appears controversial, lead- 
ing to the questions of why the r exp -divergence vanishes 
in a portion of gapped systems despite well-defined en- 
ergy gap formation at the Fermi level and whether the 
bottleneck concept is completely invalid in these gapped 
systems. Theoretical studies have been unable to satis- 
factorily answer these two questions over the last decade. 



In this Letter, we develop a theory of photo-excited 
carrier relaxation dynamics with the objective of resolv- 
ing the abovementioned issues. We postulate that the 
lack of T oxp -divergence even in gapped systems can be as- 
cribed to the presence of transverse acoustic (TA) phonon 
modes. Further, we state that an ensemble of TA modes 
in gapped systems serves as a high-capacity thermal sink, 
and energy release to the sink from other phonon modes 
facilitate efficient cooling of carriers, in other words, a 
significant reduction in r exp close to T c . The proposed 
theory is validated by the quantitative agreement with 
experimental data of r exp for peanut-shaped Ceo poly- 
mers, a typical charge density wave (CDW) material that 
does not show any carrier relaxation anomalies. 

First, we briefly review the conventional bottleneck 
concept for r oxp -divergence. It is based on the as- 
sumption that photo-excited carrier relaxation in gapped 
systems is regulated by the anharmonic slow decay of 
phonons [l|, [2(| . By the absorption of a pump laser pho- 
ton, electrons in the valence band are excited far above 
the initial states, and they rapidly accumulate in the up- 
per end of the gap through carrier-carrier and carrier- 
phonon interactions. Then, the carriers emit high-energy 
phonons (HEPs) whose energies are higher than the gap 
width 2A to relax into the lower end of the gap. The 
HEPs produced can be reabsorbed to create new carriers 
above the gap [see Fig. [TJa)] or they can decay in an an- 
harmonic manner into low-energy phonons (LEPs) that 
no longer excite carriers because their energy is lower 
than 2A [see left-hand side panel in Fig. Hfb)]. If the re- 
absorption probability per unit time is much larger than 
the inverse of the anharmonic decay time (t _1 ), photo- 
excited carriers and HEPs settle in nearly steady states 
that obey the Fermi and Bose distribution functions, re- 
spectively, with temperature T' that is higher than the 
lattice temperature T [l[. Meanwhile, LEPs obey the 
Bose distribution function with T, because they remain 
unperturbed after the laser pulse incident. Consequently, 
the carrier relaxation is dominated by the energy transfer 
from HEPs to LEPs, that is, r oxp ~ r, which is described 
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FIG. 1: (Color online) (a) Emission-reabsorption process of 
HEPs across the energy gap with width 2A. (b) Diagram of 
three-phonon scattering processes through which HEPs dissi- 
pate. Only the scatterings that involve TA modes are relevant 
to the lack of divergence in r. (c) Phonon density of states of 
TA modes with cutoff frequency Q,ta, and those of LA modes 
with Slh . The maximum frequency of the LEP is Hi (see text 
for its definition). 



plies that TA modes serve as a thermal receiver into 
which HEPs can dissipate, as a result of which the r- 
divergence vanishes even at T c . Below, we prove that 
this holds true in certain gapped systems. 

To formulate the anharmonic decay time r of HEPs, 
we consider the time evolution of the phonon distribu- 
tion function. Following the above discussion, we divide 
phonon excitations into two groups. One group consists 
of HEPs (tii < ojq.LA < ^h) in equilibrium temperature 
at T' that can participate in the reabsorption process. 
The other group involves TA modes (w qi TA < ^ta) and 
LEPs [wq,LA < 2A/H(= Cli)] in equilibrium at T that 
do not participate in the reabsorption process. Here, 
Wq.TA and Wq LA represent the phonon dispersion rela- 
tions, where q is the wavevector. tif^tii, and Ota are 
the cutoff frequencies for the HEP, LEP, and TA modes 
within the Debye approximation, respectively The dis- 
tribution functions for the two groups are given by 



n(u q ,j) = 



exp 
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(1) 



with an appropriate variable x = 

The rate of change in n(uj q j) due to phonon-phonon 
collisions is written as 21] 



by the time evolution of the two temperatures, that is, 
T = T'(t) and T = T(t). 

It should be emphasized that in the conventional bot- 
tleneck concept, only longitudinal acoustic (LA) phonon 
modes are considered. Here, we point out the unnoticed 
but important role of TA phonon modes in the HEP's 
decay [see right-hand side panel in Fig. [TJb)]. Because 
TA phonon modes generate no density modulation in the 
lattice, they cannot interact with photo-excited carriers 
within the deformation potential theory. This fact im- 
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dt 



= J col [n(uj q j)]. 



(2) 



The collision integral J C o\[ n { UJ q ,j)] describes three- 
phonon scattering, and it is defined by 



J co i[n(wo)] = 



where 
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H^/l 2 (^Sa + Sb) , (3) 
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Sa = {[n(u)o) + l]n(wt)n(u)2) — n(wo)[n(u)i) + l][n(w2) + l]}5(fou)o — tvui — fru2), 
Sb = {[n(uJo) + l][n(ui) + l]n(w 2 ) — n(uj )n(uJi)[n(uJ2) + 1]} 5(fiw 2 - fi^o - hui). 
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Here, oj s (s = 0,1,2) is an abbreviation of C0q s ,j s , and w Qi ta < ^ta, respectively. The energy Ej is defined by 
wi^f is the matrix element between the initial and the 
final state [22j ]. 

Note that J co \ regulates the time evolution of the total 
energy through the relationship 

d 



— (£ , LA + -Bta)= ^ ^qXA^col + ^2 ^q.TAJcol, 



Ej = ^fiw q jn(uj q j) = / hojn(u>)pj(uj)duj, (7) 

1(3) J ° 

where y = £!ta(^;) for j =TA (LA), pj is the density of 
states defined by 
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where the first and second summations on the right-hand where z = J1ta(^/i) for j =TA (LA), vy is the num- 
side run over qs satisfying < cj q .LA < tii and < ber of phonon branches participating in the three-phonon 
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scattering, N is the number of unit cells, and 9 is the 
Heaviside step function; the total density of states is 
given by p — pta + Pla [see Fig. [He)], satisfying the 
normalization condition 1 = J p(uj)du/N . The Debye 
approximation we have used in Eq. ([5]) is valid when 
k B T < M7ta(< K^h), which holds for many gapped sys- 
tems below T c . We simplify n{uS) in Eq. ([7]) as ksT/hoj 
[H to obtain Ej — CjT, where Cj = fcBOj2/ 3 /3 (ks is the 
Boltzmann constant). In addition, we focus on the fact 
that each sum in Eq. ([6]) gives a non-zero value only when 
one or two HEPs contribute to the three-phonon scatter- 



ing represented by J co \ [see Fig. 
a result, Eq. ([6]) is rewritten as 



ljb) as an example]. As 
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using the definitions 
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/la - w 2 ll0 {V A T + V B r), (11) 
Ita = w 2 2l i(V c T + VdT') + wIj 2 VeT, (12) 



where 7fc = 2 7 rfc|a^ /c a^ A /(S 3 7V 2 ). T is given by 
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where e is the dimensionless photoexcitation energy [lj. 
To derive Eqs. (fTTj) and ([12]), we assumed that the matrix 
element is momentum independent, that is, w^t = wi 
(or W2) =const. when TA modes are absent from (join 
in) the three-phonon scattering. Vx (X = A to E) in 
Eqs. (ITU and (fl"2"]) are functions of SI^^ta, and Q/j. 

Equation (fTTJ|) is the main finding of this study, because 
it clearly shows the contribution of TA phonon modes to 
carrier relaxation. Figure Ufa) shows the T-dependence 
of r that we have formulated in Eq. ([TO"]). The parame- 
ter p (= W2/1V1) is the coupling strength through which 
the HEPs decay into TA modes, and it is tuned from 
to 0.5 in increments of 0.1. With an increase in p, the 
magnitude of r decreases over the entire T range. When 
p is smaller (larger) than ~ 0.3, r increases (decreases) 
as T approaches T c from below. The most significant 
phenomenon is the drastic reduction in r at T c with in- 
creasing p. In fact, r at T c is inversely proportional to 
p 2 , as shown in Fig. [U[b). These results indicate that 
the lack of r-divergence at T c is attributable to efficient 
phonon-phonon coupling between the HEP and the TA 
mode characterized by p. 

To explain the microscopic mechanism for the lack of 
divergence, in Fig. HJc), we show the T-dependence of 
the magnitude of /la and /ta given in Eqs. (UT]) and 
(fT2]) . Here, /la(ta) quantifies the efficiency of the HEP's 
energy dissipation through the anharmonic interaction 
with LA (TA) phonon modes. In the limit of T ->• T c 
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FIG. 2: (Color online) (a) Numerical result of r(T) based 
on Eq. ()10[) . See text for detailed numerical conditions, (b) 
Inverse square law of r at T = T c with respect to p. (c) 
T-dependence of /la and Ita defined by Eqs. and (|12l) . 
Because Tta differs from zero at T = T, the r-divergence at 
T is strongly suppressed. 



(that is, A —5- 0), /la vanishes (irrespective of p) but /ta 
converges to a finite value as long as p ^ 0. Therefore, 
we obtain a non-diverging r at T c if p ^ 0, which readily 
follows from Eq. (fT0| . When p — 0, on the other hand, 
/ta = for arbitrary T [because /ta oc w 2 ; see Eq. ([T2]) ]. 
In this case, /la = /ta = at T c so that r — > 00 at T c . 
We thus conclude that the anharmonic decay of HEPs to 
TA modes plays a prominent role in the efficient cooling 
of HEPs as well as in determining the significance of the 
T-divergence at T c . 

Upon reducing the temperature to zero, the energy 
dissipation of HEPs gradually reduces due to monotonic 
decreases in /la and /ta [see Fig. M,c)]. The decrease in 
/la(ta) is attributed to the reduced phonon population, 
which results in a monotonic increase in r at low T, as 
shown in Fig. dja). A similar increase has been found 
in various gapped systems 0-0, lOp 12, 15, 18 1, and it is 
attributable to the inefficient cooling of HEPs. 

Now, we apply the proposed theory to photo-excited 
carrier relaxation in quasi one-dimensional Cgo poly- 
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mers. It was previously observed in experiments 
that the Ceo polymers undergo the CDW transition 
at T c = 60K, forming a well-defined energy gap at the 
Fermi level. This result implied the possibility of r oxp - 
divergence at 60K; nevertheless, optical pump-probe in- 
vestigations [Hj of the Ceo polymers revealed a mono- 
tonic variation in r cxp near T c , whose origin has yet to 
be clarified. This problem is solved by considering that 
the twisting phonon modes of the Cgo polymer [261 ] play 
the same role as the above-described TA modes. Fig- 
ure [3] shows the numerical reproduction (indicated by 
lines) of the experimental data 18j (indicated by circles) 



of the T-dependent T CX p of the Cgo polymers. An over- 
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FIG. 3: (Color online) Experimental data of r cxp in the Ceo 
polymer (after Ref. [18(]) and their numerical reproduction 
based on Eq. JEjJ. The data for T < T c = 60K fits the 
theoretical curve with p — 0.3. 



all agreement between the theory and the experiments 
is obtained by assuming p ~ 0.3. The other parame- 
ters we used are fixA = 220cm -1 , £l h — 360cm -1 , and 
2A(0) = 360K; the first two values were estimated from 
the phonon model for ID Cqq polymers (26[, and the last 
value gives 2A(0) /ksTc = 6 consistent with many CDW 
compounds }27| . 

The generality of Eq. (fT0|) is of great importance. It 
is applicable to r cxp in other gapped systems such as 
Tl-based cuprate superconductors 0, Hf and the solid 
fullerenes K3C60 and RbsCgo [3j- Even these materials 



show a lack of r oxp -divergence; however, no theoretical 
studies have attempted to clarify their relaxation anoma- 
lies. We believe that the proposed theory will serve as 
a unified framework for nonequilibrium carrier dynamics 
in gapped systems. 

In conclusion, we have developed a theory to describe 
the energy dissipation from HEPs (fiw q .LA > 2A) to 
LEPs (ftWg,LA < 2A) and TA modes in gapped systems 
below T c . This theory enables the evaluation of the In- 
dependence of r as a function of the coupling strength 
p = W2/W1 between the HEPs and the TA modes, and 
it explains the crossover between the diverging and the 
non-diverging behaviors of photo-excited carrier relax- 
ation. The latter behavior can quantitatively account for 
the anomalous carrier relaxation time in Cgo polymers, as 
measured by recent pump-probe laser experiments. The 
coupling between the HEPs and the TA mode also sug- 
gests the variation of r exp -divergence in typical gapped 
systems. 
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